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Abstract 

The decays 1:^,1], "q' 77 are investigated up to next-to-next-to leading order in the 
framework of the combined \/Nc and chiral expansions. Without mixing of the pseu- 
doscalar mesons the Nc independence of the tt*^ and r/ decay amphtudes is shown to 
persist at the one-loop level, although the contribution of the Wess-Zumino-Witten term 
to the pertinent vertices is not canceled by the Nc dependent part of a Goldstone-Wilczek 
term. The decay amplitude of the singlet field, on the other hand, depends strongly on 
N(. and yields under the inclusion of mixing also a strong Nc dependence for the ry decay. 
Both the rj and i]' decay are suited to confirm the number of colors to be Nc = 3. 
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1 Introduction 



Experimental evidence suggests that we live in a world with three colors. At high energies the 
Drell ratio for e~^e~ annihilation supports = 3, while in the low-energy regime below 1 GeV 
the anomalous decay of the 7r° into two photons , tt^ —>■ '-ff, is presented as a textbook example 
to confirm the number of colors, see e.g. |Tj. The quark charges are assumed to be independent 
of Nc yielding at tree level a width Tj^o^^^ proportional to N^, thus being quite sensitive to 
the number of colors. However, it was shown recently in |2] that the cancellation of triangle 
anomalies in the standard model with an arbitrary number of colors leads to dependent 
values of the quark charges. For three light flavors {u, d, s) Nc enters as a quantized prefactor 
of the Wess-Zumino-Witten (WZW) term 011], but the vertex with one pion and two photons 
is completely canceled by the Nc dependent part of a Goldstone-Wilczek term jSJEl. A similar 
cancellation also occurs for the decay i] 77, if one neglects 77-77' mixing. 

On the other hand, the quark triangle diagram of the microscopic theory describing the 
decay r/o —>■ 77 of the flavor singlet is Nc dependent and hence due to rj-r]' mixing the decay 
width of the 1] will also pick up an Nc dependent portion. If one works with different up- and 
down-quark masses, rriu 7^ md, the 7r° will also undergo mixing with the rj-rj' system and its 
decay width into two photons will have an — albeit small — Nc dependent piece. 

Loop diagrams which have not been discussed in [2] can be another source of A^^c-dependence 
for the two-photon decays. As we will see, the vertices of the one-loop diagrams contain indeed 
an Nc dependent piece that does not cancel out in the sum of the WZW and Goldstone-Wilczek 
terms. 

In order to investigate systematically the effects of mixing and loops for the two-photon 
decays of 71^,1^ and rj', we include the t]' explicitly within the combined framework of chiral 
perturbation theory (ChPT) and the l/Nc expansion, so-called large Nc ChPT jHllZllH]-^ In 
this theory, the r/o is combined with the octet of pseudoscalar mesons (vr, i^', T^g), since in the 
large Nc limit the axial U{1) anomaly vanishes and the r]' converts into a Goldstone boson. In 
the present work, we evaluate the decay amplitudes of 7r°, 77 and r]' up to next-to- next-to-leading 
order at which loops start contributing in large Nc ChPT. 

The paper is organized as follows. In the next section, we discuss the WZW term under the 
inclusion of the t]'. We will see that it can be decomposed into the conventional SU{3) WZW 
Lagrangian, the Goldstone-Wilczek term and counter terms of unnatural parity. In Sec. 3 the 
calculation for the decays up to next-to-next-to-leading order is presented. Numerical results 
and the importance of Nc dependent contributions are discussed in Sec. 4. Sec. 5 contains 
our conclusions and the scaling behavior of the coupling constants under changes of the QCD 
running scale is presented in App. El 

2 Wess-Zumino-Witten term 

In this section we will first briefly outline the method of extending the SU{3)Ft x SU{3)l 
chiral rotations of the effective Lagrangian in conventional ChPT to ^^(3)^ x U{3)l in a more 
generalized framework including the r/' (HI ISj- Within this approach the topological charge 

^For an alternative approach to include the rj' without employing large Nc counting rules, see, e.g., [HI- 
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operator coupled to an external field 6 is added to the QCD Lagrangian 

c = c^ao - -^^e{x)ii,{G,Mn (1) 

with G^y the gluonic field strength tensor, G^y = e^yaisG""^ its dual counterpart, and trc is 
the trace over the color indices. Under f/(l)/j x U{1)l the axial U{1) anomaly adds a term 
— -^;^N f{aR — aL)tic{Gfj,yG^^) to the QCD Lagrangian, with Nj being the number of different 
quark flavors and ur/l the angle of the U{1)r/l rotation. The vacuum angle 6{x) is in this 
context treated as an external pseudoscalar source that transforms under axial U{1) rotations 
as 

9{x) e'{x) = e{x) +i\ndeiR- In det L (2) 

with R,Le U{1), so that the term in the anomaly proportional to the topological charge 
operator of the gluons is compensated by the shift in the 6 fleld. There are, however, further 
axial anomalies which are accounted for within the effective theory by the WZW term. The 
dynamical variables of the effective theory are the pseudoscalar mesons (vr, i^, r^g, r^o) that live 
in the coset space U{3)r x U{3)l/U{3)v = U{3). They are most conveniently collected in a 
unitary matrix U{x) G U{3) with a phase given by 

detf/(a;) = e^^(^-). (3) 

The fleld describes the singlet fleld rjo and is the extension from the standard framework 
where the effective fleld is an element of SU{3). Under chiral rotations, the effective fleld U{x) 
transfroms as 

U'{x) = R{x)U{x)L\x), (4) 

so that its phase changes by 

^{j\x) =iIj{x) -i\n det R + \n det L. (5) 

Hence, the combination ip = ip + 6 remains invariant in the effective theory under chiral 
U{3)r X U{3)l rotations. The covariant derivative of U involves left- and right-handed sources 

D^U = d^U - ir^U + lUl^i (6) 

with r^ = + and = v^ — a^. In terms of these building blocks the WZW effective action 
is given by jHl ll| 

+iU dl I U^r + idrUl U^r - U^rUl + T^^U^dr U I - S^f/tr U I 
+J:J dl + T.^dl I - + is J E^/ - -{R^ L), (7) 

where = WdU and we adopted the differential form notation of |6j, 

V = dx^Vfj_, a = dx^a^, r = v + a, I = v — a, d = dx^d^j,. (8) 

with the Grassmann variables dx^ which yield the volume element dx^dx'^dx°'dx^ = e'^'^'^^d^x. 
The brackets (...) denote the trace in flavor space and the operation {R ^ L) indicates the 
interchange of r with / as well as of U with W, so that, e.g., is replaced by = UdU'^ . 
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In order to extract the SU{3) version of the WZW term, it is convenient to introduce the 
notation 

U = e^^U, detU = e~'^. (9) 

As the field ip = ip + 6 is gauge invariant, U transforms in the same manner as U under chiral 
rotations and its covariant derivative is defined as 

D^tJ = d^U - i{v^ + a^)U + iU{v^ - a^), 

af, = ttf,- l{a^) - ^d^e = a^- lD^,9. (10) 

In in] it has been shown that the WZW term can be decomposed as 

Swzw{U,v,a) = S^vzw{U,v,a) + j B (11) 

with 

B = ^(ijUFrDUDU^ + iFjDWDU + 2FrUFiW + 2F^ + 2F-^) 

1447r2 ^ ' r / / 

+ l^{Fr- - Fi) {Fr - Fj) - zDOiFrDUU^ - FjU^ DU)) (12) 

and 

Ff = df- if'^ , Fi = dl - iF . (13) 

The quantities f, I are the QCD renormalization group invariant parts of the left- and right- 
handed gauge fields r = f + ^D9, 1 = 1 — \D6 with DO = d9 + 2{a). The left- and right-hand 
side of Eq. (fTT|) actually differ by two contact terms which transform in a nontrivial manner 
both under chiral rotations and under the QCD renormalization group. In order to obtain a 
renormalization group invariant anomaly, one must remove these contact terms ,.6j. Since these 
two terms involve the singlet axial vector field (a^) and the derivative of the QCD vacuum 
angle, df^6, they are not relevant for the present work and can safely be neglected. 
The first term in Eq. (jlll) contains the WZW term for the SU{3) effective theory 

j d^x Cwzw{U,v,a) = Swzw{U,v,d), (14) 

while the second one is gauge invariant and does not contribute to the anomaly. It is straight- 
forward to show that the expression B can be absorbed by contact terms of unnatural parity 
at fourth chiral order 

d^x Cpi = ill i> {FrDUDU^ + FiDU^DU) + 2L2 {F^UFiU^) 
+2L3 %l) {F^ + Ff) + lUDO^FrDUU^ - FiU^DU) 
+2L5iJ [{Fr){Fr) + {Fi){Fi)'^+2Le^ {Fr){Fi), (15) 

where we employed the notation 

F^. = dr- ir^, Fi = dl - it^. (16) 
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The vacuum angle 6 has served its purpose and will be omitted for the rest of this section. 
Since we are interested in radiative decays, we will furthermore set the external vector and 
axial- vector fields 

r = l = v = -eQA (17) 

with A being the photon field. The anomalous Lagrangian C^^zw in Eq. ()14|) relevant for the 
two-photon decays at the one loop level reduces then to 

j d^x {tr^UUvUv + t^vdv + t^dvv -itlv) - {R^ L) (18) 



487r2 

with U = e^^U and = U^dU. 

The quark charge matrix Q of the u- d- and s-quarks has usually been assumed to be 
independent of the number of colors with Q = ■|diag(2, —1, —1). However, the cancellation of 
triangle anomalies requires Q to depend on j2] 

Q = diag(Q„.Q,.0.) = idiag(i- + 1. -1 - 1. i- - l) = Q + (l - f (19) 

with Q = idiag(2, — 1, — 1) being the conventional charge matrix, while the second term is 
proportional to the baryon number and gives rise to the Goldstone-Wilczek term. The anoma- 
lous Lagrangian of Eq. (fT^ decomposes into the conventional WZW Lagrangian of the SU{3) 
theory with the charge matrix Q and a Goldstone-Wilczek term which vanishes for Nc = 3 



SwzwiU, v) = SwzwiU, v) + (l - (20) 



with V = —eQA and 
SwzwiU , v) 



487r2 
iNce- 



487r^ 



j {i±l-±l)Q)A 

-J {2i±, -±^)Q'' + Q i±,U^QU -±^UQU^))dA A, (21) 



SawiU,v) = ^^ J {n)A-^, J {i±,-±,)Q))dAA. (22) 

It has been shown in that the Nc dependent part of the Goldstone-Wilczek term cancels 
both the 7r-27 and the 77-27 vertices of the WZW Lagrangian, yielding at tree level a decay 
width for these decays which does not depend on A^^^, if one neglects rj-r]' mixing. However, at 
the one-loop level other vertices involving kaons will contribute to the decays. One can easily 
show that, e.g., the vertex with two photons and vr"*", vr", tt*^ of the WZW term is canceled by 
the Nc dependent piece of the Goldstone-Wilczek term, in agreement with the observation that 
the number of colors does not appear in the effective theory for two flavors [2]. The vertices 
involving kaons, on the other hand, do not cancel and an A^^ dependent piece remains for the 
vertices. Consider as an example the vertex with two photons and vr", K~^, K~ that contributes 
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to the decay 7r° — 77. The WZW term yields the vertex (neglecting mixing of the vr" with the 
rj-rj' system) 

-^^K+K-dn^dAA, (23) 
whereas the Nc dependent piece of the Goldstone-Wilczek term leads to 



367r2 



K+K-dn^dAA. (24) 



Clearly, both terms do not compensate and a dependence on A^^ remains in the final expression 
for the vertex. It is therefore of interest to study the N^. dependence of the two-photon decays 
at the one-loop level. 



3 Radiative decays at one-loop order 

In the framework of large Nc ChPT the expansion in powers of momenta and light quark masses 
is combined with the l/N^ expansion by ordering the series according to 



P 



OiVS), m, = Oi6), l/Nc = Oi6). (25) 



In this bookkeeping, the WZW term Swzw is of order 0{6), whereas the one-loop diagrams of 
the decays involve the ratio rriq/ p with / ~ 0{V6) being the pseudoscalar decay constant in 
the chiral limit and are thus of order 0{6^), i.e. next-to-next-to-leading order. 
Our starting point is the WZW effective action of the U{3) theory 



SwzwiU, = J d X C^yzwiU, v) 

= -^^ / d^x {T.^U^dvUv + ^^vdv + T.^dvv- iT.lv) -{R^L) (26) 
We expand the quark charge matrix Q in powers of 1/Nc 

= 0'°' + <3'", (27) 

where the superscript denotes the order in the combined large and chiral counting scheme, 
i.e. Q^°^ (Q^^^) is of order 0{1) {0{S)). With U = exp(z0), one obtains for the three decays 
from S-^yzw the tree level contributions 

d'x = -^{d<P Q') AdA = -^(0 [(Q^°^)' + 2gWg(i) + (g(i))2]> dA dA (28) 

since Q^'^\Q^^^ commute with the diagonal entries of (p — 0^,0^,0''. The terms on the right- 
hand side of Eq. contribute at orders 6, 6^ and 6^, respectively, if one disregards the 
A^^c-dependence of in which a factor 1// = 0{l/\/I^) has been absorbed. Hence, within large 
Nc ChPT the 0^ and (p^ decay amplitudes start at order 0(5^), whereas the decay amplitude 
for the singlet field 0° is of order 6. 
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At fourth chiral order the unnatural parity Lagrangian consists of more terms, which are 
gauge invariant, see Eq. (fTH|) . 

with 

+ {{F,){Fr) + {Fi){Fi))+2V,{^){Fr){Fi) + ..., (30) 

where we have presented only the contact terms which contribute to the decays at the order 
we are working. 

The potentials Vi are odd functions in t/S, so that the leading contribution in the l/N^ 
expansion is linear in ip. The Lagrangian CpA can be expanded in powers of 1/Nc 

44 = 4^)+ 4? + ..., (31) 

where the superscript denotes the order in the 5 expansion with 

d^x 4? = 2L2 {FrUFiU^) + 2L3 ^ {F^ + Fl) (32) 

and 



d^x = 2L5 [{Fr){Fr) + {Fi){Fi)') + 2U ^ {Fr){Fi). (33) 

! from £'•4'* and C^' 
J) 

-eQ dA, we obtain for £*■ 



The contributions from C^"^^ and C^^4 are of order O(p^) and 0{^p'^), respectively. Setting 



d'x C^> = 2e'[L, + 2L3]V^((g(°))^ + 2Q^'^Q^'^ + (Q^'yf) dA dA 
2^ /3 1 



e%{--^)^dAdA (34) 



where ki = 2{L2 + 2L3) and the last term proportional to (Q*-^-*)^ has been omitted, since it is 
of order (9(5^) and thus beyond our working precision. In a similar way, the terms from £^^4 
reduce to 

d^x 4? = e2fc2^(Q(°)) (Q(°)) dA dA (35) 

with k2 = 2{2L^ + Lq). 

From the renormalization group invariance of the effective Lagrangian it follows that ki and 
/c2 transform as (c/. App. [XJfor details) 

A;^^" = ZAk2, (36) 

where Za is the multiplicative renormalization constant of the singlet axial current = 
1^7^ 75g which transforms as (A^)*"^" = ZaA^^ under changes in the QCD running scale. 
At sixth chiral order the relevant terms for the decays read 

d^'x 4 = zWi{ij){Ux^F^ + x^UFf) + iW2{i^){x^F,UFi + Ux^UFiU^F^) 
+ m{^P){Ux^){F^ + Fl) + iW,{'iP){Ux^){U^FrUFi) 

+ iW^{i^){Fr + Fi)i\F,.U + UFi\x^) + h.c. . (37) 
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The quark mass matrix Ai = diag(mu, rud, rris) enters in the combination x = 'iBM. with B = 
— (0|gg|0)//^ being the order parameter of the spontaneous symmetry violation. Expanding in 
powers of l/N^ one obtains 

= Cf + (38) 

with the contributions C^^^ at 0{NcP^) 

d'x Cf> = twf\[Ux^ - xU^]F^ + [x^U - U^x]F?) 

+ iwf ([x^ - U^xU^]FrUFi + [Ux^U - x]FiU^Fr) (39) 

and order 0{p^) 

d'x Cf = ) ^ {[Ux^ + xU^F^ + [x^U + U^x]F?) 

- w'i^ ^{[x^ + U^xU^FrUFi + [Ux^U + x]FiU^Fr) 

+ iwf\Ux^ - xU^){F^ + Fl) + iwf\Ux^ - xU^){U^FrUFi) 

+ iwf\Fr + Fi){FrUx^ - U^FrX + U F^^ - FiU^x). (40) 

where the potentials Wi have been expanded according to Wi = wf^ + iwl^^tp + 0{tjj'^). 
The explicitly symmetry breaking terms reduce to the structures 

d^x = A;3e2(0x[(Q^°^)' + 2Q^^^Q^'^]) dA dA (41) 
with = —4:{'wf^ + W2^) and 

d'x £f = e'(k,^x{Q^''>r) + h{<Px){{Q^''>r) + h{Q^''>){<PxQ^''>)) dA dA (42) 



with k^ = —A{w^^ + W2^), k^ = —2{2wf^ + wf^), k^ = —8wf \ respectively. The scaling law 
for the parameter k^ is given by (c/. App. \^ 

kr = ZAk, + ^[ZA-l]k,, (43) 

while the remaining parameters k^, k^ and k^ remain put. 

Having discussed the tree diagram contributions to the decays, we now turn to the calcula- 
tion of the loops at order 5^. After expanding the WZW Lagrangian in the meson fields the 
contributing pieces at one-loop order read 

2 



d'x C^,^ = [0,Q2]]) _ (d0[0,Q][0,Q]))A dA 

iN e 

-^{d(Pd(Pd(PQ)A + ... . (44) 

The mesons inside the loops do not undergo mixing, as 0^, cj)^ and (jP loops do not contribute. 
The first two terms in Eq. ()44|1 contribute via tadpoles, whereas the last one represents a 
vertex of the unitarity correction, corresponding to Figure lb. As the diagonal components of 
commute with Q, tadpoles with 0^, 0^ and 0° do not contribute, and since the photon couples 
only to charged mesons, the unitarity corrections arise due to charged meson loops. 
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(a) (b) 
Figure 1: One-loop diagrams contributing to — 77. In (b) the crossed diagram is not shown. 

At order 0{6^) only the Q^^'' piece in Eq. (jlij) contributes. From = \l it follows then 

that the first term vanishes and the tadpoles are entirely due to the second term. Performing 
the loop integration, one obtains for the decays the tadpole contributions 



(45) 



i=3,8,0 

with the tadpole given in dimensional regularization by 

d'^l 



{2ttYP -ml + ie 
where L contains the pole at d = 4, 



2L + 



167r2 



In 



(46) 



/"V 1 1 



1671^ \d- 4: 2 
and /i is the regularization scale. The coefficients d^ read 



(ln47r + r(l) + l)) 



(47) 



d^ 



0, 



V2' 



dt 



d 



K 



3' 
1 



dl = - 



d° 



2 

7!' 
2 



(48) 



Evaluating the unitarity corrections at order 0{5^) for on-shell photons and replacing Q in 
Eq. (gH) by Q^, since the contribution from Q*^^) is beyond the order we are working, one 
obtains exactly the same contribution as for the tadpoles but with opposite sign. Hence, the 
one-loop corrections to the decays at order 0{5^) compensate each other and the first non- 
vanishing non-analytic piece will show up at order 0{5'^). This is also in agreement with 
previous calculations in conventional ChPT, in which the chiral logarithms were compensated 
completely by wavefunction renormalization and replacing / by the physical decay constant 
in the tree level expression jTUlEI- However, within large N^. ChPT the cjr' and 0^ decay 
amplitudes start at order 0(5^), so that the leading non-analytic corrections to the physical 
decay constant and wavefunction renormalization will contribute at order (9(5^) and do not 
affect the amplitude up to order 0{5^). Any divergences from the loop diagrams discussed 
above could then only be renormalized by counter terms of the Lagrangian of unnatural 
parity. This would be clearly in contradiction to previous results IT^ . 
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So far, the (j)^ amplitude does not contain an explicit Nc dependence due to the WZW term, 
but both the 0^ and 0" fields undergo mixing which results in the mass eigenstates 7] and rj' . 
Here, we work in the isospin limit of equal up- and down-quark masses, m = niu = ma, so that 
the (p^ field decouples from the 77-77' system. In the next section, we will give an estimate on 
isospin breaking effects due to different quark masses niu ^ nid- 

In order to describe 77-77' mixing up to one-loop order, one must take into account the 
following terms of the effective Lagrangian 

= + (49) 

which reads at lowest order 

C^^) ^ t{D^U^D>^U) + ^{xU^ + Ux^) - \Ti,\ (50) 

The next-to-leading order Lagrangian C'^^^ — 0{5) contains the terms 

£(1) = L,{D^U^D>^U{x^U + U^x)) + L^{x^Ux^U + U^U^x) 

+LA,D^^D^i; + i^A2^{x^U - U^x), (51) 

and is given by 

£(2) = L,{D^U''D^^U){x''U + C/tx) + L^ix^'U + U'^x? + Mx^U - U'^x? 

+iL,sD^7j;{D''U''x - D'^Ux'') + iL^^i^i^^xUh - x'Ux'^U) + 0{N,p^). (52) 

The terms r]Q{D^W D'^U) and 7^0 (x^^ + ^X^) have been omitted in C^'^\ since the pertinent 
unknown coupling constants represent OZI violating corrections. Moreover, as indicated in 
the last equation, counterterms of order 0{NcP^) with new unknown coupling constants will 
contribute at order 5^. We will neglect these contributions throughout, assuming that they are 
of small size and do not alter our conclusions. 

The fields (p^ and 0° are related to the mass eigenstates rj and rj' via 

08 = [cos7?(i) - sin7?(0U(^)] V + ^ [sin^?^'^ + cost^^^U^^^] 77' 

(fP = [cos 7?(^U(') - sin 7?(2)i3] ^ + [sin d^^) A^^^ + cos d'^^^B] 77'. (53) 

The decay constants and F^, are defined by 

(0|w5A°g|77')=i^/2p^F° (54) 
with the normalization (A" A'') — S^^, while the angles 7?^*^ correspond to the mixing angle up 
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to order 0{6^) which arises in the diagonahzation of the 0^-0° mass matrix 
sin 279(°) : 



4a/2 — 
3 m^, - m2 



sm 



2^9(^) = sin2t9(0) ( 1±Al + 1^[2lP _ )](m2, - m^) - ^L^r 



16 

n 

-|i [m^ - m\\ {2Lf + SL^) r (1 - ^A,) + [m^ - m^] [llm^ + ml] \ 

+ ^ K - <] [5< + m^] (24^(1 - 1a,) - lSJ) - ^ m>« + -1 m?,Ai^) 
^ AM (m^ [| cos^ + ^ sin 2^(°) + 2 sin^ ^(o)] 



.ft' 



[4 cos^ ^(°) + 2V2 sin 279(°) + 2 sin^ ^(°)]) 



^ A J) (m^ [| sin2 ^(0) - -1= sin 2^(o) + 2 cos^ ^W] 



3^2 ^""^12 " V2 ' 

TT 

-m^[4sin2^(°) - 2\/2sin2^(°) + 2cos2^(°)]) + 0{N^p^)^ (55) 

where ^ 93 MeV is the pion decay constant defined in a similar way as in Eq. (jSU, 
and m^,m^/ are the diagonal entries of the 77-77' mass matrix. Furthermore, h'f''^ and A^'^) = 
m^/(167r2) ln(7n|/ jj?) are the finite parts of the LECs and the loops, respectively, after renormal- 
ization. It is straightforward to verify that the angles t?'-*-' do not depend on the regularization 
scale of the effective theory. 

The quantities 771,^ and ttik denote the pion and kaon masses at leading order 

ml = 2Bm = ml{l - ^ml[2Lt^ - L^] - A(2< + ml)[2Lt^ - L^] - ^A« 

+ 4;^[cos2 ^(0) + 2 sin^ ^W] AW + 4;^[sin2 ^(0) ^ 3 cos^ d^'''^]^}^)] , 



ml = B{m + m,) = ml[l - j-,ml[2Lt^ - 4^] - ^(27^^ + ml)[2Lt^ - L^i 



^ cos27?(o)aW) + ^— [4771^ cos^ 79(0) - (37n2, +m2)sin2 79(°)]A5'^ (56) 



3F2 J 12F2 
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The expressions A^^^ and B read 



a(2) 4^2 2 



2Lr^ + 3lS 



+ 



;i + Ai)i/^ 



2L^'^ + 3L 



(r) 
18 



(r) 



-(r) 
18 



Ai 



64 



9Fi 



(57) 



For the 0^ decay we have only kept the pieces up to next-to-leading order, since the terms 
beyond that order contribute at 0{5^). In the case of the 0° decay, on the other hand, one 
must keep also the contributions at next-to-next-to-leading order. The values of the couplings 
Ai, A2, and L25 are not known and depend on the running scale of QCD. As they represent 
OZI violating corrections, we will omit them, but, strictly speaking, we cannot expect that all 
neglected terms vanish at the same scale. Furthermore, the parameter r is related to the mass 
of the 7]' in the chiral limit which was estimated in to be about 850 MeV. This translates 
into a value of r ^ 1 x 10~^ GeV^. 

Including the mixing from Eq. (|55|l we obtain the amplitudes 



|l7rO + ( |^[cos^(i) - sin^9W^(i)] ^ ^^^^^ ^(1) j^(2) _ gin^(2)^] 



V 



+ (|^[sin^(i) + cos^9(°U(^)] + ^[sin^9(^U(2) + cos^?^^)^] 



where tt" is related to the 0^ field via 

with the pion Z-factor 
and the decay constant 



dA dA, 



(58) 



/ 



4 



™2 r(r) 



(59) 

(60) 
(61) 
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to the order we are working. The coefficients B^p read 



Bt, 



:{XsQ') + V2h{Xsx[iQ^'^r + 2Q^'^Q^'^]) 



:{XsQ') + V2h{Xsx[iQ^'^r + 2Q^'^Q^'^]) 



4V27r2 

+V2h{Xsx){{Q^''>r) + V2h{Q^'^){XsxQ^''>) 



Brji 



+Jlh{x){{Q^'r) + J\h{Q^'^){xQ^'). 



(62) 



Due to the N^. dependence of the quark charge matrix Q the expressions B^^ and B^ start at 
order (5^, whereas B^,' contains a piece of order 0{5). Substituting these relations into Eq. (j58|) 
yields the decay widths 



r 



r„/ 



||[cost9(^) - sin^(°U(^)] + ^[cos^^'U^'^ - sin^^^^B] 
+ cos 19^°^^^^^] + ^[sini?^^^^^^^ + cos-i^^^^iS] 



(63) 



with a = e'^/Air. In the 6 expansion the leading order contribution to the decay width of the rj 
is given due to mixing by the leading contribution in 5^/, and is comparable in size with the 
Bjj portion. The numerical values will be discussed in detail in the next section. 



4 Numerical analysis 

Equations ()55p and (j37|) are utilized to obtain values for the mixing angles ■i?*^*-' and the expres- 
sions A^^\ B, respectively. For the LECs L^^l g g, Ly, we take values which follow from matching 
the U{3) theory to the SU (3) framework by integrating out the singlet field at the regular- 

ization scale /i = 1 GeV, L['\fj.) = -0.5, 4''^(/i) = 1.0, 4''^/^) = -0.3, Lj = -0.3, L^s\l^) = 0-7 
(all in units of 10~^). Note that integrating out the singlet field only alters the LECs L^q\Lj 
and Lg'^^ and their values in the U{3) framework are within the phenomenologically determined 
error ranges of the SU{3) LECs. Moderate variations of these LECs yield small changes in 
the decay amplitudes with the largest changes induced by variations in Lj roughly at the 
10% level. Our conclusions are therefore not altered, if slightly different values for the LECs 
are employed. Using the experimental values for the pseudoscalar meson masses, we obtain 
■(9(0) = -21.8°, = -15.8°, ^ _i9.8°. 
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By employing Eq. (jHS)), we can now fit the ratios Bj,/ Fj,, Br,/ F^, Bn' / F^, to the decay widths 
FttO^^^, r^^^^, Vr^i^-y-f. The experimental values for the decay widths are |20j 

r^o^^^ = 7.74 ± 0.55 eV, 
r^^^^ = 0.465 ± 0.045 keV, 

r^,^^^ = 4.28 ± 0.34 keV, (64) 

and the fit to the central values yields 

B^/F^ = -0.133 GeV-\ 
Brj/F^ = -0.0522 GeV-\ 

= -0.192 GeV-^ (65) 

For the pion decay constant F-j^ we employ the physical value F^^ ^ 93 MeV, while F^ can be 
extracted from a one-loop calculation with F^ = 1.34F^ It is consistent to take the one-loop 
results for and F^, since the difference with respect to the next-to- leading order expressions 
shows up at 6^ in the decay amplitude and is, therefore, beyond our working precision. The 
values for B^/F^^ and B^^/F^ from the fit are close to the contributions from the anomalous 
WZW term, B^^'^/F^ = -0.136 GeV-\ B^^'^/F^ = -0.0587 GeV-\ indicating that the 
portions from the counter terms of unnatural parity are small. Omitting higher orders beyond 
S^, they contribute with a relative strength of about 2% to B-j^ and 10% to S^. In order to get 
an estimate for F^,, it is thus justified to assume that the counter term combination in 5^/ is 
small as well. Of course, both the counter terms and F^, depend on the renormalization scale 
f^QCD, but we will assume that for a certain range of figco the counter term contributions are 
negligible. For such /Iqcd, the ratio B^,^^ /F^, is then reproduced by setting ^ 1.16 F-^, a 
value slightly larger than in previous calculations [TIH lllj. 

In particular, we would like to investigate, whether a clear statement can be given on the 
number of colors by utilizing the l/iV^ expansions of the decay amplitudes. The cancellation of 
Witten's global SU{2)l anomaly requires to be odd |22j. The standard model with Nc = l 
is without strong interactions. We will therefore restrict ourselves to a comparison of the 
numerical results for Nc = "i and Nc= 'b. Setting all non- anomalous contact terms of unnatural 
parity to zero, the decay width for the in a world with A^^ = 5 reads F^^=^ = 1.002 keV to 
be compared with the decay width in the real world with three colors, F^'^i^ = 0.511 keV. For 
the T}' we obtain F^;^^^ = 4.21 keV and F^^^^ = 12.8 keV. The experimental values for the t] 
and r}' decays clearly rule out A^'c = 5 and varying the values for the omitted coupling constants 
of the counter terms within realistic ranges does not alter this conclusion. 

Finally, we would like to give an estimate on the dependence of the 7r° decay width due 
to different up- and down-quark masses. In the case of different up- and down-quark masses, 
the 0^ field undergoes mixing with both the 0^ and 0° field. In order to get an estimate, we 
will restrict ourselves to the mixing at leading order in the 5 expansion. The fields 0'^, 0^ and 
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6° are then related to the mass eigenstates via 
6^ = ^fr:^- 



^ ( cos {r] + evrO) + sin ^^^^ + e'7r°- 

V 

V2 



- sin {r] + e7r°) + cos ^9^°) (//' + eV") ) (66) 



with the mixing parameters 



\/3md-mu 

5 



eo- 



rris — m 
cos^9(°) - y2sin^9(o) 
1 + ^tan?9(o) 



sin^^W + y2cos?9(°) 
' = l-;i3Ctg^(0) • 

The parameter eo can be expressed in terms of physical meson masses by applying Dashen's 
theorem which implies the identity of the pion and kaon electromagnetic mass shifts up 
to C(eV) 

There have been estimates in the literature that Dashen's theorem is significantly violated 
at higher orders due to chiral symmetry breaking effects jEl ^1 1^ • On the other hand, a 
recent non-perturbative approach to the hadronic decays of r] and rf indicated that higher order 
corrections to this low-energy theorem may be small ^Hl- In any case, the isospin- violating 
effects in the decay widths constitute a small correction, so that it is safe to employ Dashen's 
theorem in the present work. 

This time the fit to the central experimental values yields 

B^jF^ = -0.134 GeV"\ 
B^/F^ = -0.0598 GeV-\ 

Br,'/F^, = -0.208 GeV"^ (69) 

and setting = 1.07F^ in the ratio B^,^^ j F^, reproduces the fitted value. 

It should be emphasized that we fitted our results to the current world average value for 
r^o^^^ = 7.74±0.55 eV [201 • On the other hand, it is possible to estimate in a model-dependent 
way the size of the counter term contributions. In [21, e.g., the values of the counter term 
contributions to the 7r° decay have been estimated by means of a QCD sum rule for the general 
three-point function involving the pseudoscalar density and two vector currents. Within that 
approach the authors find a slightly enhanced width of F^o^^^ = 8.10 ± 0.08 eV. 

A comparison of the numerical results for A^^c = 3 and A^c = 5 and setting all non-anomalous 
contact terms of unnatural parity to zero yields F^''^^ = 8.00 eV, F^''^^ = 8.18 eV for 
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the pion decay, T^,^!^^ = 0.456 keV, T^,"'^^ = 0.920 keV for the r] and T^Zjj = 13.7 keV, 
r^!^^^ = 4.28 keV for the rj'. Again, A^c = 5 is ruled out by the rj and rj' decays, but no rigorous 
statement can be made for the vr^, although the value for A^^c = 3 is in better agreement with 
the current world average. 

5 Conclusions 

In the present work, we investigated the two-photon decays of vr", t] and t]' in the combined 1/Ac 
and chiral expansions. The cancellation of triangle anomalies in the standard model requires 
the quark charges to depend on N^. We have shown that the WZW term of the U{3) effective 
theory decomposes into the conventional anomalous SU{3) WZW Lagrangian, a Goldstone- 
Wilczek term and counter terms of unnatural parity which involve the singlet field rjo. The 
independence of the 7r° and rj decay amplitudes on Ac which was shown in [2^ to occur at tree- 
level due to partial cancellations of the WZW term with a Goldstone-Wilczek term, persists at 
one-loop order, although the vertices of the pertinent loop graphs do exhibit an dependence. 

We performed a one-loop calculation including counter terms up to next-to-next-to leading 
order in large ChPT in which 77-77' mixing has also been taken into account up to one-loop 
order. Within the bookkeeping of large Nc ChPT, the leading contribution to the 77 decay arises 
from mixing with the 77'. 

From a fit to the experimental decay widths and under the assumption that higher orders 
beyond our working precision can be neglected it follows that contributions from the counter 
terms are small. Since the cancellation of Witten's global SU{2)l anomaly requires A"c to 
be odd and a world with Nc = 1 has no strong interactions, we compare the cases Ac = 3 
and A^c = 5. The numerical results of the rj and 77' decay widths for A^c = 3 are close to 
the experimental values and clearly rule out the case Ac = 5. We have furthermore given an 
estimate on the A'c dependence of the pion decay due to different up- and down-quark masses 
by taking Tr^-rj-i]' mixing at leading order into account. The A^c dependence of the vr^ decay is 
smaller than the experimental uncertainty and is therefore not suited to extract the number 
of colors. We conclude that both the 77 and the 77' decay show clear evidence that we live in a 
world with three colors. 

It has been pointed out in |2j that at tree level the process 77 — > 7r+7r^7 is proportional to 
A"^ and should replace the textbook process 7r° — > 77 lending support to Ac = 3. It will be 
thus of interest to investigate the decays 77 7r^7r~7 and 77' tt^ti^^ within the framework 
of large A^, ChPT 
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A Scaling behavior of the coupling constants 

In this appendix, we derive the scaling behavior of the coupling constants which contribute to 
the decays. The transformation properties of the constants L2,3 have already been given in jSj 
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and we merely quote the result here. The quantities L2 and L3 are renormalized according to 

Lr = ZaU - It = ZaU - (A.l) 

with 

.= ^M^. (A,2) 

1447^2 ^ > 

Since the singlet field ^ scales as ^^^"^ = Z^^^p, the scaling behavior of Lg*^^ ensures that the 
Lagrangian Cy^zw + ^^p4 remains invariant to order 5^ under changes of the QCD running scale. 

~ ~(3) 

In order to study the transformation properties of L5 g we rewrite £^4 

= 2^ [2L5 + Lq] {dv) (dv) + 2tp [2L5 - Lg] {da) {da) (A.3) 
This yields the transformation properties 

{2L, + L,)''^ = Za{2L, + L,), 

(2Z5 - L,)- = Z\{2L, - U) - ^^^3^~^\ (A.4) 

so that the Lagrangian C^zw + ^^^p} + '^^4'' remains renormalization group invariant. 

On the other hand, the contact terms Wi and W2 in the Lagrangian of sixth chiral order 
can be written as 

iWi{^)e^'^{Ux^F^ + x^UFf) + iW2{^)ei^ {x^ FrU Fi + Ux^UFjU^Fr) + h.c. + .... (A.5) 



The ellipsis in Eq. ()A.5|) denotes terms with more than one fiavor trace which involve contri- 
butions from other contact terms and are irrelevant for the discussion of the scaling behavior 
of Wi and W2. From (jA.5j) we obtain the transformation properties 

ty2(x)"'^" = W2{ZAx)e-^^^^-^^\ (A.6) 

Expanding the potentials Wi in the singlet field ip yields for the two leading expansion coeffi- 
cients 

= ZAwf^ + ^(^^ - ^ = 1' 2. (A.7) 
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